Mon. Not. R. Astron. Soc. 000,[T]-?? (2012) Printed 10 January 2013 (MN I*TbX style file v2.2) 



The merger rates and sizes of galaxies across the peak epoch of star 
formation from the HiZELS survey. 

2 John P. Stott 1 *, David Sobral 2 , Ian Smail 1 , Richard Bower 1 , Philip N. Best 3 , 
O James E. Geach 4 

Institute for Computational Cosmology, Durham University, South Road, Durham, DH1 3LE, UK 
(■H i 2 Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands 
^ , 3 SUPA, Institute for Astronomy, Royal Observatory of Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK 

■ 4 Department of Physics, McGill University, Ernest Rutherford Building, 3600 Rue University, Montreal, Quebec, H3A 2T8, Canada 



O 

u 

Oh! 

O 
& 

<N 
> 

m 
oo 

<N 

i— H 

<N 



X 



ABSTRACT 

We use the HiZELS narrow-band Ha survey in combination with CANDELS, UKIDSS 
and WIRDS near-infrared imaging, to investigate the morphologies, merger rates and sizes of 
a sample of Ha emitting galaxies in the redshift range z = 0.40 — 2.23, an epoch encom- 
passing the rise to the peak of the star formation rate density. Merger rates are estimated from 
space- and ground-based imaging using the M20 coefficient. To account for the increase in 
the specific star-formation rate (sSFR) of the star forming 'main- sequence' with redshift, we 
normalise the star-formation rate of galaxies at each epoch to the typical value derived from 
the Ha luminosity function. Once this trend in sSFR is removed we see no evidence for an 
increase in the number density of star-forming galaxies or the merger rate with redshift. We 
thus conclude that neither is the main driver of the enhanced star-formation rate density at 
z ~ 1 — 2, with secular processes such as instabilities within efficiently fuelled, gas-rich discs 
or multiple minor mergers the most likely alternatives. However, we find that ~ 40 — 50% 
of starburst galaxies, those with enhanced specific star formation at their epoch, are major 
mergers and this fraction is redshift independent. Finally, we find the surprising result that the 
typical size of a star-forming galaxy of a given mass does not evolve across the redshift range 
considered, suggesting a universal size-mass relation. Taken in combination, these results in- 
dicate a star-forming galaxy population that is statistically similar in physical size, merger rate 
and mass over the ~ 6 Gyr covered in this study, despite the increase in typical sSFR. 
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1 INTRODUCTION 

The peak in the volume averaged star formati on rate for galax- 
ies occurs in the redshift range z — 1 - 3 dLillv et al. 1996; 
iMadau et alJl 19961 : 1 Sobral et alj|2013h . At this epoch, the star for- 
mation rate (SFR) in typical gala xies is an order of m agnitude 
higher than in the local Universe ( Reddv & SteidellE009h . This is 
the era when most of the stars in the Universe were formed and rep- 
resents the peak in black hole activity. The task is now to address 
'how' and 'why' the Universe was so different then. 

A picture is emerging in which the dominant mode of star for- 
mation at this earlier epoch is very different to that in the local Uni- 
verse. Rather than the quiescent formation of stars that is the norm 
in today's Universe, violent episodes of star form ation are domi- 
nated by the formation of super-star clusters (e.g. ISwinbank et alj 



* E-mail: j.p.stott@durham.ac.uk 



2010b). However, the origin of these differences is somewhat con- 
troversial: one picture, which has some observational support, is 
that they are driven by an increase in the galaxy merger rate (e.g . 
Somerville eTai1l200li : iHopkins et aTll2006l : IConselice etaD 2003. 
2008), but other theories have suggested that it is the result of the 
higher rate of gas accretion expected in the high-redshift Universe 
(iKeres et al.ll2005l : iDekel et alJl2009l) . It is therefore important to 
study the SFR, merger fractions and gas content of these galaxies 
in order to identify the processes responsible for driving this epoch 
of enhanced activity. 

In recent years the presenc e of a star forming 'main- sequence' 
seen in the local Universe (e.g. iBrinchm ann et al. 2004) ha s been 
confirmed at increasing ly high redshift dElbaz et al.l I2007L EoTH : 
iDaddi et alJ2007l : lRodighiero et al.l201ll : ISargent et alJ2012h . This 
is a relation between SFR and stellar mass for star forming galax- 
ies, with a typical specific star formation rate (sSFR, the ratio of 
the star formation rate to the stellar mass of the galaxy) found to 
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increase with redshift (lElbaz et al.ll201lb . Galaxies that lie off this 
relation with sSFRs too high to be in the typical star-forming popu- 
lation are often described as 'starbursts' and a re thought to be trig- 
gered by violent events such as major mer gers (IHopkins et al. 1 2006: 
lElbaz et al.ll2007l : lRodighiero et alfcOllh . 

From a theoretical perspective, in the A Cold Dark Matter 
(ACDM) paradigm dark matter halos merge hierarchically from 
the bottom up, with the largest halos created at later tim es (e.g. 
lLacev & Coiell 19931 : ICoTe et al.ll200d : ISpringel et al.ll2005l) . As the 
galaxies trace the underlying dark matter we therefore expect those 
to merge hierarchically also. However, it has been known for some- 
time that the most massive galaxies appear to hav e older stel- 
lar populations than thei r less massive counte rparts (iCowie et alj 
1 19961 : Bower et alJl200d iGilbank et ai1l2010h . Environment also 
plays a key role with massive quiescent galaxies typically liv- 
ing i n denser en vironments than lower mass star-forming galaxies 
toresslej[l98Q ). There are several ways to reconcile these obser- 
vations with hierarchical merging which are implemented in phe- 
nomenological, semi-analytic models that seek to reproduce obser- 
vations of galaxy evolution by populating dar k matter halos from 
N-body simulation s with mock galaxies (e.g. iBower e t al. 2006; 
ICroton et aDl2006h . A reasonable match is achieved through in- 
teractions and feedback mechanisms that cease star formation in 
massive galaxies within massive dark matter halos, requiring that 
these galaxies build up their stellar mass at late times by so called 
'dry' mergers which trigger no significa nt new star format ion due 
to the lack of available cold gas toe Lucia & BlaizotH2007l) . 

In the high-redshift Universe the cold gas fraction in galaxies 
is higher tha n at low-redshift and thus there is more fuel for star for- 
mation (e.g. lTacconi et al.ll20ld : lGeach et al.ll201ll) . It is therefore 
possible to more easily tri gger significant star-forming events dur- 
ing mergers (iSomerville et al.ll200ll) or through hig h gas accretion 
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Cacciato et al.l2 012). The latter process leads to the intriguing pos- 
sibility of the enhanced star-formation rates at high redshift being 
dominated by secular evolution rather than mergers. In fact while 
some obs ervations suggest an increase in the merger fraction with 
redshift dConselice et al .1 12003) others seem to prefer in-situ galac- 
tic processes ov er galaxy-galaxy merging, or at lea st a mixture of 
these processes dLotz et al l200a lElbaz etaLll2007b 

To test whether it is galaxy mergers or secular processes that 
dominate and drive galaxy evolution at the peak era for star for- 
mation, a method to distinguish between galaxy mergers and non- 
mergers needs to be implemented. The two main methods of es- 
timating the merger fraction are counting close pairs of galax- 
ies, under the assumption that they will subsequently m erge (e.g. 
iLe Fevreetai1l2000l : iLin et al.l [20081 : iBluck et alJl2009h . and us- 
ing a method of ide ntifying galaxies with a merging morphology 
(e.g.lConselice et alJl2003l: lLotz et aljl2004 IConselice et all 2008 ; 
lLotz et al.l2008l : IConselice et alJl2009h . The results of both of these 
methods often suggest that the merger fraction increases with red- 
shift and, depending on the mass range considered, the merger 
fraction at z > 1, where the star formation rate density peaks, is 
roughly ~ 0.1 — 0.3 on average (but with some systematic offsets 
between studies) compared to a fraction < 0.1 in the local Uni- 
verse. A third, potentially more reliable, method is to employ de- 
tailed integrated field unit observations of z — 1—2 galaxies to look 
for merger signatures in the dynamics of the galaxies. Such stud- 
ies, although generally smaller in sample size, also find a merger 
fracti on of ~ 0.3 (e.g. lForster Schreiber et al.ll20Q9t Tshapiro et all 
120081) . 



In order to study the star-forming population, an excellent 
tracer of ongoing star formation is the Ha emission line which 
is less affected by dust obscuration than shorter wavelength star- 
formation tracers (e.g. UV continuum light or [OII]3727). Un- 
fortunately beyond z = 0.4, Ha is redshifted out of the opti- 
cal window, thus high redshift studies of star formation have been 
limited to either using the obscuration-effected short wavelength 
tracers or studying small samples of Ha emitters using conven- 
tional near-infrared spectrographs. However, in the last few years 
panoramic narrow-band surveys have started to provide large sam- 
ples of Ha- selected gala xies (e.g. the High-red s hift (Z) Emission 
Line Survey HiZELS. iGeach et al .1 |2008l 120121 : iGarn et all 1201 Ol: 
Sobral et a l. 2009, 2o7b U2012ll2013l and the studies of IVillar et alJ 
2008 and lLv et aDl201 lh . Narrow-band surveys provide a well un- 



derstood, volume- selected sample of star-forming galaxies allow- 
ing for straight-forward analysis of trends with SFR, mass and size 
etc. They provide emission line information over large areas of the 
sky and are thus able to probe a significant range of the Ha lumi- 
nosity and stellar mass functions for star-forming galaxies, required 
for a n unbiased analysis of the star formation r ate density (SFRD, 
e.g. lGeach et aDl2008l ; [Sotoal et al.ll2009U201 1120131) . This selec- 
tion method has also been shown to be extremely effective at de- 
tecting intrinsically faint galaxies, helping to overcome the bias to- 
wards massive galaxies associated with photometric redshift selec- 
tion. 

In thi s study we use the z — 0.4 — 2.23 HiZELS sample pre- 
sented in I Sobral et al.l d2013l) . to not only analyse the merger rate 
as a function of redshift and stellar mass but also as a function of 
the well-determined SFR. We can therefore test whether it is ma- 
jor mergers that drive the rise to enhanced activity seen at these 
epochs. In contrast to earlier studies, which analyse Hubble Space 
Telescope (HST) rest frame UV morphologies, with the advent of 
the WFC3 camera we can also study the rest-frame optical bands 
for a subsample of our galaxies that lie within the CANDELS re- 
gion of our survey and use this to calibrate morphologies derived 
from deep, wide-field, ground based near-infrared imaging, better 
matched to the extent of the full HiZELS fields. We also analyse the 
size-mass relation for star-forming galaxies over this epoch in or- 
der to study the size evolution which may also indicate the merger 
history of these systems. 

The structure of this paper is as follows. In 32 we describe the 
HiZELS narrow band sample and the imaging data. We then derive 
SFR for the sample and analyse the evolution of the number den- 
sity of galaxies above a given SFR. The size-mass relation is then 
studied in order to look for an evolution. A method for automating 
morphological classification is defined and this is used to study the 
merger rates of the galaxies in our sample and how they evolve and 
depend on SFR and mass. Finally, we discuss our findings in the 
context of understanding the physical processes that occur within 
galaxies, that lead to the rapid downturn in the global volume aver- 
aged SFR below z ~ 1. 

A ACDM cosmology (ft m 0.27, Q A 0.73, H 70 
kms -1 Mpc -1 ) is used throughout this work and all magnitudes 
are AB. 



2 THE SAMPLE AND DATA 
2.1 The HiZELS survey 

HiZELS dGeach et all 120081 : ISobral et all 120131) is a Campaign 
Project using the Wide Field CAMera (WFCAM, ICasali et al.1 
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2007) on the United Kingdom Infra-Red Telescope (UKIRT) and 
exploits specially designed narrow-band filters in the J and H 
bands (NBJ and NBH), along with the H 2 S1 filter in the K band, 
to undertake panoramic, moderate depth surveys for line emitters. 
HiZELS targets the Ha emission line redshifted into the near- 
infrared at z = 0.84, 1.47 and 2.23 using these filters. In addition, 
the UKIRT data are complemented by deeper narrow band obser- 
vations with Subaru Suprime-Cam NB921 imaging dSobral et alj 
l2012Ll2013h to obtain Ha emitting galaxies at z = 0.4 and the 
[Oil] emission from the z = 1.47 Ha sample, as well as deeper 
WFCAM and Very Large Telescope near-infrared imaging through 
the H2SI filter in selected fields. The survey is designed to trace 
star-formation activity across the likely peak of SFR density and 
provide detailed information about a well-de fined statistical s am- 
ple of star-forming galaxies at each epoch (see lBest et al.ll201ob . 

In this study we concentrate on the main HiZELS sample 
of z = 0.4, 0.84, 1.47 and 2.23 Ha emitters in both the UKIRT 
Infrared Deep Sky S urvey, Ultra Deep Survey (UKIDSS UDS, 
lLawrence et all2007l Alm aini et al. in prep.) an d The Cosmic Evo- 
lu tion Survey (COSM OS, IScoville et aDl2Q07h fields as described 
in Sobral et al 1 d2013h and we refer the reader to that paper for full 
details of the catalogues used. These data cover areas of 0.6 — 1.6 
square degrees depending on the field and waveband. The narrow 
band excess sources are visually inspected to remove image arte- 
facts and, to ensure the galaxies are at the desired redshift, spectral 
energy distribution (SED) fitting and optimised colour-colour se- 
lections are used to provide clean sam ples of Ha emitters in the 
four redshift slices (Sobral et al . 2013). The excess narrow -band 
flux is then converted into an emission line luminosity. For the 
analyses in this paper we take these cleaned catalogues and in- 
troduce cuts to ensure that the data in each narrow-band filter are 
complete to the same flux limit across the entire area observed. 
These final catalogues contain: 428 Ha emitters at z = 0.40, 595 
at z = 0.84, 420 at z = 1.47 and 372 at z = 2.23 down to the SFR 
limits ~ 0.2, 3.0, 12.0 and 25.0 M©yr _1 respectively (assuming 
Aua = 1.0), to an Ha equivalent width lower limit of 25A. 

The star formation rates for the HiZELS sampl e are cal- 
culate d from the Ha luminosity and the relation of iKennicuttl 
dl998h (SFR(M©yr _1 ) = 7.9 x 10~ 42 L( H a) (ergs" 1 )) 
assum ing a dust extinction Ah a — lmag (see ISobral et all 
2013). Stellar masses are computed by fitting SEDs 
to the rest-frame UV, optical and near-infrared data 
available (FUV, NUV, U,B,g,V, R, z, 7, z, Y, J , H,K 

3.6/im, 4.5/im, 5.8/xm, 8.0/xm coll ated in [S obral et al. I2013L 
see references therein), following ISobral et al.l (1201 lb and the 
reader is referred to that paper for more deta i ls. The SED templates 
are generated with the Bruzual & Charlotl d2003h package using 
Chariot & Bruzual (2007, unpublished) models, a lChabrierl(l2003h 
IMF, and an exponentially declining star formation history with 
the form e _t//r , with r in the range 0.1 Gyrs to 10 Gyrs. The 
SEDs were generated for a logarithmic grid of 200 ages (from 0.1 
Myr to the maximum age at each redshift being s tudied). Dust 
extinc tion was applied to the templates using the ICalzetti et al.l 
d2000h law with E(B-V) in the range to 0.5 (in steps of 0.05), 
roughly corresponding to Ahq ^0 — 2. The models are generated 
wi th different metallic ities, including solar; the reader is referred 
to ISobral et al] d201 lh for further details. For each source, the 
stellar mass is computed as the median of stellar masses of the 1 a 
best-fits over the range of parameters. 

In Figure \T\(left) we plot the number density of galaxies, from 
the combined UDS and COSMOS fields, above a stellar mass of 
10 10 M Q and a given SFR, against redshift. From this we can see 



that for a given SFR limit the number density increases rapidly 
with redshift. This is a manifestation of the fact that a typical star- 
forming galaxy has a greater sSFR at higher redshift, forming stars 
more rapidly for a given mass. In order to look for trends with 
redshift we therefore define a quantity which we term the epoch- 
normalised star formation rate ENSFR which is the SFR of a galaxy 
divided by the SFR* (z). SFR* (z) is the star-formation rate derived 
from the quantity L^ a found by fitting a Schechter function to the 
Ha luminosity fun ction at a given redshift, which we take from 
ISobral et al.l d2013h . We note that normalising the SFR to SFR* (z) 
accounts, to first order, for the increase in sSFR with redshift. How- 
ever, significant evolution in either the slope of the SFR - stellar 
mass relation or the dust obscuration would invalidate this. 

The values of SFR* essentially double for 
each HiZELS redshift interval considered with 
SFR* - 7.0, 14.0,29.0, and 57.0 Moyr -1 for z = 
0.4, 0.84, 1.47 and 2.23 respectively. Interestingly, this same 
behaviour is seen in the evolution of the typical sSFR from 
lElbaz et all bOllh with sSFR - 0.2,0.4,0.8, and 2.0 yr -1 , 
again at these redshifts. We suggest that this is because the Ha 
luminosity (and thus SFR) function evolves significantly more 
than the stellar mass function. 

In Figure Q] (right) we plot the number density of galaxies of 
a given mass above the thresholds SFR/SFR*0) = 0.6, 1.2, 2.4. 
From this plot one can clearly see that the number of star-forming 
galaxies with their SFR normalised to the typical SFR at that epoch 
is broadly constant. This means that the number density of star- 
forming galaxies of a given mass and ENSFR does not evolve sig- 
nificantly over the period studied here. This demonstrates that the 
star-forming population is constant with redshi ft but simply evolv es 
in sSFR. This is similar to the result found in ISobraletal.ld2013l) in 
which there is no strong evolution in the Schechter parameterisa- 
tion of the normalisation of the Ha luminosity function, 0H a . We 
discuss the implications of this in §5\ 



2.2 Imaging data 

In this study we analyse near-infrared imaging from the 
space-based HST/WFC3 Cosmic Assembly Near-infrared Deep 
Extragalactic Legacy Survey (CANDELS, iGrogin et all 1201 lb 
iKoekemoer etaDboilh and the gro und-based UKIDSS UDS and 
the WIRCam Deep Survey (WIRDS. lBielbv et all2012h . 

The CANDELS imaging we use is from WFC3 F160W cover- 
ing a 2-orbit depth over 720 sq. arcmin of the UDS. The CANDELS 
imaging has a pixel scale of 0.06 arcsec and a point spread func- 
tion (PSF) with a FWHM of 0.18". The CANDELS data are well 
suited to this project for which we require high resolution imaging 
in the rest-frame optical, however to obtain the wider area cover- 
age needed to build up a statistical sample of rarer high-mass sys- 
tems from HiZELS we also need to use ground-based near-infrared 
imaging. 

The UKIDSS UDS K-bmd imaging covers an area of 0.8 
square degrees, to a depth of K = 24.6 (5cr, AB) with a pixel 
scale of 0.13 arcsec and a PSF FWHM of 0.7". The WIRDS K- 
band imaging covers a total effective area of 2. 1 square degrees and 
reaches an AB 50% completeness limits of ~ 24.5 across the COS- 
MOS field, it has a pixel scale of 0.15 arcsec and a PSF FWHM of 
0.7" and is thus comparable to the UKIDSS UDS. 

The combination of these three near-infrared imaging datasets 
allows us to probe the rest-frame optical morphologies and sizes 
of the HiZELS galaxies over a wide range in luminosity while at 
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Figure 1. Left: The number density of HiZELS galaxies above a stellar mass of 10 10 M© and a given SFR plotted against redshift. The SFR > 25 MQyr -1 
lines are offset slightly in z for clarity. Right: The number density of> 1O 1O M0 galaxies above an epoch normalised star formation (ENSFR) threshold. 
We define ENSFR as the ratio of SFR to SFR* (z) (with SFR* (z)derived from the Lj^ a ie the typical SFR from the Ha luminosity function at that redshift, 
ISobral et al.l2013h . In this way we remove the trend that the average sSFR of galaxies increases with redshift. As there is no evidence of a significant trend this 
demonstrates that the number density of typical star-forming galaxies does not evolve significantly with redshift and thus the increase in the SFRD is purely 
an effect of increased typical sSFR. 



the same time providing a rest-frame optical view of the galaxies' 
stellar distribution. 



3 ANALYSIS 
3.1 Sizes 

Before studying the morphologies and the merger rates of the 
galaxies in the HiZELS sample, we first assess their typical sizes. 
This is interesting from a galaxy evolution perspective, as an in- 
crease in size with cosmic time may imply that mass is being built 
up either through mergers or accretion or that the mass is being re- 
distributed somehow. If there is no direct evolutionary connection 
between the galaxy populations at each epoch then changes in typ- 
ical size may suggest differing formation scenarios. Importantly, it 
will also help us to understand the reliability of the morphological 
classification as the smallest galaxies will be most affected by the 
resolution of our ground-based imaging. 

The surfa ce photome try of galaxies is often described by a 
Sersic profile (SersiclQi 



I(r) - 7 e exp < -b n 



l/n 



(1) 



where J(r) is the intensity, r is the radius from the centre of 
the galaxy, r e is the scale radius, I e is the intensity at r e , n in the 
exponent is a free parameter widely known as the Sersic index and 
b n = 2n — 0.327; a coefficient chosen so that r e is the half-light 
radius defin ed as the radius wh ich encircles half the light from the 
galaxy (e.g. Graham etal.lll996h . 

To measure the sizes of the galaxies we fit a 2-dimensional 
Sersic profile to th e galaxy images using the GALFIT (version 3) 
software package dPeng et al.ll2002h . This software requires rea- 
sonable initial input parameters such as position, apparent magni- 
tude and ellipticity, all of which are estimated b y first running the 
SEXTRACTOR package teertin & Arnouts|[l996l) so that the itera- 
tive fitting process converges to the correct solution in the shortest 



possible time. GALFIT deconvolves the point spread function which 
is dominated either by the telescope itself, in the case of HST, or 
by the atmospheric seeing for the ground-based imaging. To this 
end we check that the effect of seeing has been correctly accounted 
for in the analysis of the ground-based imaging by comparing the 
CANDELS derived sizes to those from the UKIDSS UDS imaging 
for the same galaxies. Figure |2] shows this comparison of galaxy 
sizes for a sample drawn from a combi nation of all fou r HiZELS 
redshift slices and a sample of BzK jDaddi et alj|2004l) galaxies 
in the UDS field (the photometry to select BzK galaxies is taken 
from the UDS catalogues, Almaini et al., in prep). These two in- 
dependent size measurements are correlated and scattered around 
the 1-to-l line with Ar e /r e ~ 0.4, which confirms that the sizes 
recovered are comparable, demonstrating that GALFIT is able to 
successfully account for the seeing. 

We note that there may be some selection effects and biases in 
size measurements, in that galaxies with large-sizes can be missed 
due to low surface brightnesses and co mpact galaxies may have 
sizes overestimated ( B ar den et al .1 120051) . The former is less likely 
as the HiZELS galaxies are selected on their Ha emission. Also, 
Figure |3 demonstrates that there is no significant bias in size esti- 
mates between the ground and space-based analysis of the smallest 
galaxies so we take this as evidence that their sizes are not overes- 
timated. 

Figure [3] shows the size-mass relations at each redshift slice. 
We perform linear fits to this relation of the form log 10 r e = 
a (log 10 (M*) — 10) + 6, where r e and M* are in units of kpc and 
M respectively and we normalise the fits to M* = 1O 1O M . Ta- 
ble Q]contains the results of these fits at the four redshift slices con- 
sidered. From these fits we find the surprising result that the typical 
size of a star-forming galaxy with log 10 M* = 10 does not evolve 
significantly out to z = 2.23, with r e = 3.6 ± 0.3 kp c on average. 
These r esults are in good agree ment with the trends of IB arden et alj 
(2005); Ichika waet al.1 J2012h who also find little evidence of an 
evolution in this relat ion or the typical siz e of star forming galax- 
ies. In a related study, Ka nwar et afl (l2008h find no evolution in the 
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Figure 2. The half-light radius measured from the UKIDSS UDS ground 
based imaging plotted against that from the HST/WFC3 CANDELS data at 
all redshifts. Solid line is the 1-to-l line. The open and filled circles repre- 
sent BzK and HiZELS galaxies respectively. The dashed and dotted lines 
represent the UKIDSS UDS PSF HWHM at z = 1.47 and z = 0.4 respec- 
tively, which bracket the other two epochs. This demonstrates that we can 
recover t he sizes of galaxi es by accounting for the ground-based PSF using 
galfit jPeng et al.l2002l) . 



Table 1. The size-mass relations at each redshift slice, of the form 
log 10 r e = a (log 10 (M*) — 10) + b. Where r e and M* are in units 
of kpc and Mq respectively. 



z 


a 


b 


r e atlog 10 (M*) = 10 








(kpc) 


0.40 


0.08±0.02 


0.55±0.03 


3.6±0.2 


0.84 


0.03±0.02 


0.54±0.01 


3.5±0.1 


1.47 


0.03±0.02 


0.59±0.01 


3.9±0.2 


2.23 


0.08±0.03 


0.51 ±0.02 


3.3±0.2 



shape of the size function of disc galaxies between 0.1 < z < 1.0 
with just an evolution in the number density of discs. However, 
other groups have found evidence for a stronger size evolution for 
the most massive (M* > 10 10 M cy)) disc-like galaxies, with a 2 — 4 
fold i ncrease in size since z ~ 2 dTruiillo et al.ll2007l : lMosleh et al.1 

limb. 

By analysing the Sersic index, n, which we obtain from the 
fitting process we divide our sample into disc-like and bulge-like 
galaxies where we define the former as having 0.5 ^ n < 2.5 and 
the latter as 2.5 ^ n < 5.0. From this we find that the fraction of 
disc-like galaxies is > 80% in each redshift slice with no evidence 
for an evolution, which is not unexpected as star-forming galaxies 
such as those s elected by HiZELS are in general found to be discs, 
consistent with lSobral et al.1 d2009h . We note that this disc fraction 
also has no trend with SFR or stellar mass. 

3.2 Morphologies 

3.2.1 Quantifying and calibrating morphology 

To quantify the morphologies of the galaxies in this study we 
choose to u se a combination of Gini and M20 coefficients first 
proposed bv lLotz et al.l d2004h . The Gini coefficient, developed by 
statistician Corrado Gini, measures the inequality among values of 
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Figure 3. The half-light radius plotted against stellar mass for the z = 0.4 
(top), z = 0.84 (upper middle), z = 1.47 (lower middle) and z = 2.23 
(lower). The solid lines are linear fits to the relations with the dotted line 
the z = 0.84 fit for reference. The dashed lines represent the PSF HWHM. 
The slope of the size-mass relation is found to be broadly constant. 



a frequency dis tribution. It wa s first applied to studies of galaxy 
morphology by I Abraham et al.l d2003l) . A Gini coefficient of zero 
expresses an equality where all values are the same (i.e. a galaxy 
with uniform surface brightness). A Gini coefficient of 1 expresses 
maximal inequality among values (i.e. where all of the flux is in one 
pixel/element). The M 2 o coefficient, describes the second-order 
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moment of the brightest 20% of pixels in the ga laxy and is sensi - 
tive to merger signatures such as multiple nuclei dLotz et al 1 2004). 
The combination of Gini and M20 can differentiate between 'nor- 
mal' star forming galaxies and Ultra Luminous Infrared Galaxies 
(ULIRGs), as well as single galaxies and merging systems. How- 
ever, there are some differences in the boundaries chos en to delin- 
eate these populations (e.g. see lLotz et al.ll200d l2008h and there- 
fore we choose to perform our own tests and calibrate the Gini and 
M20 coefficients by visual inspection. 

The Gini and M20 coefficients are calculated using 
the Gini and M20 compo nents of the galVSM software 
(iHuertas-Companv et al . 2008). This software requires a segmen- 
tation map which tells galVSM which pixels are associated with 
the galaxy. We first cutout 10" postage stamp images, taken from 
the CANDELS mosaic, around e ach galaxy and generate a segmen- 
tation map using S EXTRACTOR dBertin & Arnoutsll 996). The Gini 
and M20 codes are then run on the postage stamps and the corre- 
sponding segmentation maps. 

The sample we choose to run the initial visual inspection cal- 
ibration analysis on is that of 167 star forming galaxies in the 
redshift range 1.4 < z < 2.5 selected via the BzK method 
(iDaddi et alj|2004h which lie within the CANDELS survey region 
in the UDS. The F160W mosaic provides high resolution rest frame 
optical imaging of these galaxies. We choose this sample over the 
HiZELS narrow-band sample as it should consist of similar star 
forming galaxies but has a higher surface density and so a larger 
sample falls within the high resolution CANDELS imaging, key to 
testing the morphological classifications. 

From visual inspection of the BzK galaxy morphologies, 
the SEXTRACTOR parameters DEBLEND_MINCONT=0.1 and DE- 
tect_minarea=5 and DETECT_THRESH=1 a are found to be re- 
laxed enough to associate clear merging components of the same 
'galaxy' with one segmentation map but still stringent enough so 
as to not produce clear false positives. We note that having a DE- 
blend_mincont set too high means that unrelated galaxies would 
be considered as mergers whereas when set to a low value separate 
features within the same galaxy separate into distinct objects and 
therefore this parameter has the most effect on the M20 coefficient 
(see Appendix lAlf or a discussion of this parameter). Setting the de- 
tection threshold to low sigma values includes fainter 'sky' pixels 
in the segmentation map and thus increases inequality, raising the 
Gini coefficient. In this way one can see that the way in which the 
segmentation map is created is the most important factor in deter- 
mining the Gini and M20 coefficients and differences between how 
this is done in different studies are the reason why we choose to 
calibrate our own definitions of mergers and non-mergers. 

Using the above method, fixing the SEXTRACTOR parameters 
to those found to give the best performance, the Gini and M20 
codes are run on the CANDELS imaging with the results for the 
BzK sample are displayed in Figure ^(upper). Also included is 
a 0.35 < z < 0.45 photometric redshift sample with a similar 
magnitude range to t he HiZELS z = 0.40 sample sourced from 
Willi ams et all b009h to demonstrate that this classification tech- 
nique is not affected by redshift. 

By visually assigning the galaxies into two categories 'merg- 
ers' and 'non-mergers' with the former classification based on evi- 
dence of merging components either creating disturbed morpholo- 
gies or very close potential mergers (on- sky separation < 2"). This 
information is included in Figure [4] with the delineation between 
mergers and non-mergers found to occur at an M20 ~ —1.5 for 
both high and low redshift regimes and thus the Gini coefficient 
does not seem to add any information. Using this method there is 



a contamination of ~ 10% non-mergers in the mergers and < 5% 
mergers in the non-mergers. The simulations performed for Ap- 
pendix [A] demonstrate that the M20 coefficient is sensitive to merg- 
ing components down to a luminos ity (mass) ratio o f ^ 1 : 10 (in 
agreement with the simulations of lLotz etal As such we 

note that our analysis throughout this paper is a measure of major 
mergers only. 

The morphology codes are then run on the same galaxies 
but using the deep if -band ground-based UKIDSS UDS so we 
can compare the two independent measurements. We expect the 
higher resolution CANDELS imaging to be a truer reflection of a 
galaxy's intrinsic morphology. We also note here that it is difficult 
to measure the morphologies of the lowest luminosity galaxies in 
our sample as they tend to be smaller (see size-mass relations in 
^3. II ) and are thus more affected by the seeing of the ground-based 
near-infrared imaging. By performing tests we find that setting DE- 
BLEND_MINCONT=0.03 ensures that the M 2 o parameter selects the 
same type of mergers in the ground-based data as that derived from 
the HST data (again see Appendix |A]). The ground-based versus 
HST Gini and M 2 o values are plotted in Figure [5] By performing 
linear fits to these relations we can calibrate the ground-based Gini 
and M20 values to those derived from HST. These fits are: 



Gini candels = 0.78 Gini uds + 0.13 (2) 
M20, CANDELS = 0.68 M 20 ,uds - 0.39 (3) 

and will now be applied to the HiZELS morphologies derived 
from the ground-based near-infrared imaging. 

One potential problem with measuring the morphologies of 
galaxies at different epochs, using the same near-infrared imaging, 
is that of morphological k correction. Galaxies look smoother at 
longer wavelengths, meaning that the lowest redshift galaxies in 
our sample may artificially appear less disturbed than those at high 
redshift. When we analyse the HST Advanced Camera for Surveys 
(ACS) F814W imaging data available in COSMOS, many of the 
galaxies are very low surface brightness and therefore it is difficult 
to assess whether the morphological classifications given by the 
M20 coefficient are reliable. However, for the galaxies in the z = 
0.4 sample with Kab < 22.5, the same classifications as those 
derived from the near-infrared CANDELS imaging are recovered 
in rsj 90% of the cases, so we conclude that our results are not 
significantly affected by this. 

An additional concern is that some disc galaxies at high 
redshift are found to conta in large star-forming clumps (e.g. 
Swinbanket al. 2010b, 2012). There may therefore be a degener- 
acy between what we classify as 'mergers' and those galaxies that 
contain a small number of large star-forming clumps. It is prac- 
tically impossible to differentiate between these two populations 
without dynamical information and thus we note with caution that 
so called 'clumpy disc' galaxies may make a up some fraction of 
our 'merger' sample, if the clumps are on scales of > 4 kpc. In 
fact when we run the sub- sample of nine HiZELS galaxies which, 
from dynamical analysi s of integrated field un it data, are all found 
to contain clumps (see ISwinbank et al.l 120121 for a description of 
this sample) all of them have M20 > — 1.5 and thus we would clas- 
sify them as 'mergers'. We note here that when visually classified 
not all of these clumpy galaxies appear as clear mergers which may 
explain the non-merger interlopers with M20 > —1.5 in Figure [4] 
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visual classifications derived from the CANDELS data indicated. 
The result of analysing the morphologies of these calibrated data is 
that we now choose to delineate the difference between mergers and 
non-mergers at M20 = —1.45 which minimises the visual contam- 
ination to 22 ± 12% non-mergers in the merger region and 15 =b 7% 
mergers in the non-merger region. We note that some of the con- 
tamination of visually classified non-mergers to the merger frac- 
tion may in fact be due to galaxies with clumpy discs (see ^3.2.11 ). 
Figure [6] displays a sub-sample of the HiZELS galaxies classified 
by the M20 parameter as mergers (left) and non-mergers (right) 
for both ground and space-based imaging, with their s EXTRACTOR 
segmentation maps over-plotted. As the Gini coefficient is found to 
add little information, when using our particular analysis methods, 
Figure |7] presents a histogram of M20 values, as measured from 
the ground-based imaging of the HiZELS population at all redshift 
slices. 



• • 
• • • 

• '# o n o o 



• 9 



-0.5 -1.0 -1.5 -2.0 -2.5 -3.0 



Figure 4. Upper: The Gini coefficient plotted against the M20 value for 
the z ~ 1.4 — 2.5 BzK population and a photometric redshift sample 
with z ~ 0.4 from the F160W CANDELS imaging data in the UDS field. 
The filled red and open blue symbols are those classified as mergers and 
non-mergers respectively by visual inspection of the CANDELS imaging 
with circles representing mergers and non-mergers for the BzK popula- 
tion and squares for the z ~ 0.4. From visual inspection M20 ~ —1.5 
appears to be an excellent delineation between mergers and non-mergers. 
This demonstrates that the for our particular analysis the key parameter for 
determining whether a galaxy has a merger-like morphology is the M20 
parameter and not the Gini coefficient. Lower. The Gini coefficient plotted 
against the M20 coefficient for HiZELS galaxies at all redshifts, as mea- 
sured from the UDS K band imaging and calibrated using equations [2 and 
[3] but with morphologies visually identified from the CANDELS F160W 
image. The filled red and open blue circles are those visually classified as 
mergers and non-mergers respectively. The vertical line at M20 = —1.45 
is the value we now choose from visual inspection to delineate the merg- 
ers and non-mergers. This demonstrates that the calibrated ground-based 
near-infrared imaging can be used to derive M20 values that differentiate 
between mergers and non-mergers. 



3.2.2 HiZELS morphologies 

The number densities of galaxies in the HiZELS samples are lower 
than the BzK morphology calibration sample used in ^3.2.11 and 
therefore do not have the same level of overlap with the CANDELS 
imaging region in the UDS. We instead run the morphology codes 
on the CANDELS, UKIDSS UDS and COSMOS WIRDS imaging 
for the HiZELS samples at each of the four redshifts. The output 
Gini and M20 values for the ground-based near-infrared imaging 
are calibrated to the CANDELS values using the fits found for the 
BzK sample in ^3.2.11 As a confirmation of the calibration of the 
ground-based morphologies to those derived from the HST data the 
UKIDSS UDS Gini and M 2 o coefficients for those that line in the 
CANDELS sub-region are plotted in Figure |4](7tfwr) but with the 



4 RESULTS 

4.1 Merger fractions 

Here we define 'merger fraction' as the number of galaxies with 
a merger-like morphology (regardless of how many galaxies actu- 
ally make up this merger) divided by the total number of galaxies 
in the redshift slice. The total fraction of mergers for the HiZELS 
galaxies in the redshift bins z = 0.40, 0.84, 1.47, 2.23 are 0.33, 
0.13, 0.18 and 0.32 respectively (see Figure [7), however these are 
not comparable as they are measured for different stellar mass and 
SFR ranges at the dif ferent redshifts. 

For comparison ISobral et al.l (E009) find a higher merger frac- 
tion of 0.28 at z — .84 using the morphological classifications of 
IScarlata et al.1 (|2007) and a visual classification that included merg- 
ers and close pairs (which explains the higher merger fraction), al- 
though this is from rest-frame B— band imaging. However, when 
we study the COSMOS HST ACS imaging used in that study we 
find that many of the galaxies appear as very low surface bright- 
ness meaning that their morphological classifications are more un- 
certain. 

In Figure [8] (left) the merger fraction is plotted against stellar 
mass, with the lowest mass galaxies progressively more likely to 
be classed as mergers with ~ 5 — 20% of the star-forming pop- 
ulation being mergers at the highest stellar masses in each of our 
redshift slices. The z = 0.4, 0.84 and 1.47 trends are all remark- 
ably similar and in agreement but there is an increase in merger 
fraction at all masses to z = 2.23. However, the HiZELS selection 
is dependent on SFR, not mass and as described in $2] the typical 
sSFR for galaxies increases with redshift and therefore we need to 
investigate these effects too. 

A fraction of 10-20% mergers is seen in the most strongly star- 
forming galaxies at each redshift (Figure [8] centre). However, due 
to the flux-limited nature of the samples and the evolution of typical 
sSFR there is little overlap between different redshifts. In this figure 
the combined SFR data for all of the HiZELS redshift bins taken at 
face value may actually hint at a trend in merger fraction with SFR 
rather than any evolution with redshift (at least out to z — 1.47). 
There is some evidence of an increase in merger rate at the same 
SFR when going from z = 1.47 to z = 2.23 but again this does 
not account for the evolution in typical sSFR. 

Combining the two results above we investigate the relative 
contribution of mergers to the range of sSFR covered by our sam- 
ple, for galaxies with ENSFR > 0.2 to which we are complete 
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Figure 5. Left: The ratio of the Gini coefficient for BzK galaxies measured from the UKIDSS UDS ground based X-band imaging to that measured from the 
CANDELS HST F160W imaging plotted against the Gini coefficient measured from the CANDELS F160W imaging. The solid line is the 1-to-l relation and 
the dashed line is a fit to the observed trend. Right: The ratio of the M20 coefficient for BzK galaxies measured from the UKIDSS UDS ground based X -band 
imaging to to that measured from the CANDELS F160W imaging plotted against the M20 coefficient measured from the CANDELS F160W imaging. The 
solid line is the 1-to-l locuss and the dashed line is a fit to the relation. From these plots we can see that it is possible to calibrate the values of Gini and M20 
derived from ground-based imaging to those from HST imaging. 
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Figure 6. Left: Postage stamp images (10" x 10") from both CANDELS HST F160W (upper ten) and UDS K (lower ten) of the same HiZELS galaxies 
classified by M20 as major mergers. Right: Postage stamp images (10" x 10") from both CANDELS F160W (upper ten) and UDS K (lower ten) of the same 
HiZELS galaxies classified by M20 value as non-mergers. It is clear from this plot that mergers are well separated from non-mergers in this morphological 
classification system and that it is possible to identify mergers from the ground-based imaging. The white outlines represent the S EXTRACTOR segmentation 
maps used for the morphological analysis. 



at all redshifts, in Figure [8] (right). From this plot it is clear that 
the galaxies with the higher sSFR at all redshifts are increasingly 
more likely to have a merger-like morphology, with those with 
the highest sSFR having a merger fraction of ~ 40 — 50%. This 
strongly suggests that starbursts are more likely to be driven by 
major mergers when compared to the rest of the star-forming popu- 
lation. This is in agreement with the far infrared selected sample of 
iKartaltepe etaD d2012h who find that major mergers have, on av- 
erage, a high sSFR compared to typical star forming galaxies. It is 
also in broad agreement with results from the mass selected sample 



of iKavirai et alj d2Q12h who also find that major mergers tend to 
have high sSFR compared to undisturbed galaxies. 

We test whether the Au a = 1.0 dust correction we univer- 
sally employ is reasonable and how it affec ts our results. This is by 
including both SED fit extinction values dSobraletalJl2013h and 
those derived from the rel ation between stellar mass and extinction 
from [Gam & Bestl (120100 . We find that using these more sophis- 
ticated estimates makes little difference for the range of masses 
we study. The value of Ahq is ~ lmag at a mass of 1O 1O M 
with this value increasing/decreasing to higher/lower mass, with 
the typical range being Ah« = 0.5 — 2 mags. In fact when this 
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Figure 7. A histogram of M20 values for the four HiZELS redshift slices. 
The vertical line at M20 = —1.45 delineates mergers from non-mergers. 



more sophisticated treatment of dust obscuration is included it acts 
to strengthen our conclusions by smoothing the relations in Figure 
[8] However, we choose to keep the extinction value at Au a = 1.0 
as this is easier to compare to other works including the main re- 
sults in ISobraletalJ (120131) and to 'epoch normalise' with the Ha 
luminosity function. 

It is unlikely that HiZELS is missing a large population of 
'typical' high redshift sta r forming galaxies with high dust ob- 
scurations, as Redd v et al ] d2Q12h demonstrate that the dust con- 
tent of typical star forming galaxies actually decreases with red- 
shift. However, the HiZELS sample may miss extreme star form- 
ing and highly obscured galaxies such as sub-mm galaxies. Sub- 
mm galaxies are found to have a spread in morphologies which is 
indistingui shable from that of t ypical star forming galaxies at high 
redshift dSwinbank et al.ll2010ah and a re relatively rare objects (1- 
2xlO~ 5 Mpc~ 3 . IWardlow et al.ll201ll) . thus their omission would 
not affect our conclusions. 



4.2 Merger rates 

To calculate the merger rates (the number o f mergers per Gpc 3 per 
Gyr) we follow the prescription outlined in lLotz et aD d201 ll) : that 
the merger rate is simply the number of mergers per Gpc 3 divided 
b y the average times cale over which the merger would be observed. 
In lLotz et"aD (l201ll) this observed merger timescale is found, from 
simulations, to be ~ 0.2 Gyr, when the Gini/M2o method is em- 
ployed. We adopt this value for consistency with that study and with 
the data from other groups recalculated and used there. We note that 
as HiZELS is a narrow-band survey the volumes covered at each 
redshift slice are w ell defined with values of ~ 1 — 7 x 10 -4 Gpc 3 
dSobral et al. 1 12013). We now also assume that there are on average 
two galaxies per merger for consistency with other studies. 

For comparison with other surveys we initially cut our sam- 
ple only on stellar mass. The merger rates for galaxies with 
M > 10 9 and 10 10 M o are plotted in Fig ure [9] (left). Also plot- 
ted are values from Conselice et al ] d2003h (M B < -19 which 
approximates M > 10 9 M^) and those with M > 10 10 M 



derive d from Gini/M2o (Lotz et al]|2008l). close pairs ( Lin et al l 
2008 ) and galaxy a symm etry from IConselice et al] 120091 and 
Lopez- Saniuan et al] J2009h . These merger rates are corrected to 



the timescales calculated by Lotz et al. ( 2011) using the galaxy evo- 
lution models of lSomerville et al. (20081). 

Figure [9] (left) shows little evolution in merger rate with red- 
shift and the results are generally in good agreement with those 
found in the studies of IConselice et alJ (120031 . l20Q9h : iLin et al] 
( 2008); lL6pez-Saniuan et al.l d20Q9b where the r edshift ranges over- 
lap. The merger rates from lLotz et all (120081) are systematically 
higher, which may be because that sample is mass selected and 
therefore includes a significant contribution from merging red 
sequence galaxies which would not have been included in the 
HiZELS sample. There could also be secondary effects due to a 
mismatch in the stellar mass calculation between the studies, a 
different way of defining mergers through the M20 parameter or 
a differential in the timescales involved, so an offset is perhaps 
not unexpected. With the exception of the z — 0.4 data point, 
which is significantly higher, the HiZELS merger rates for galax- 
ies with M > 1O 9 M are also in good agreemen t with those of 
the o nly study with this approximate mass limit (IConselice et al] 
120031) . From this plot there is no strong evidence for an increase in 
the merger rate for mass-selected samples out to z ~ 2. However, 
this comparison does not account for the SFR l imits of the d ifferent 
surveys or the increase in typical sSFR with z felbaz et aljl201ll) . 

The advantage of HiZELS over these earlier studies is that it is 
unbiased with respect to stellar mass and we derive the stellar mass 
and SFR from independent measurements, i.e. SED fitting and Ha 
flux. We can therefore consider both SFR and stellar mass inde- 
pendently to split the population into sub- samples based on these 
properties. As defined in 32 we account for the increase in the typ- 
ical sSFR with redshift by employing the ENSFR. Figure ^\(right) 
shows the population split into three ENSFR bins > 0.6, 1.2 and 
2.4, for which the HiZELS observations are complete at all red- 
shifts. The first obvious thing to notice is that undulating shape of 
the plot with just a mass cut (Figure |9] left) has disappeared. In- 
stead the trends are flat, showing no evidence for an increase in 
the merger rate with increasing redshift for all masses and ENSFR 
cuts. This mass and ENSFR selected sample is a cleaner sample 
than those in Figure ^\(left) and so we suggest that the peak in the 
merger rate at z ~ 1 seen for some comparison samples may be 
due to the mixing of a mass limit with an SFR selection function 
which strongly effects photometrically- selected galaxies. 

From this merger analysis we can determine the total number 
of major mergers (with mass ratio > 1 : 10) a galaxy of a given 
mass will underg o during the epoc h covered by our study. Using 
Equation 11 from lConselicd J2006) we find that one would expect 
~ 3 mergers per star-forming galaxy with M ~ 1O 1O M0 between 
z = 2.23 and z = 0.4, or a merger every 2 Gyrs on average. We 
note that these numbers depend on the value of r the timescale over 
which mergers can be observed using the M20 method (which we 
assume to be 0.2Gvr. lLotz et alJboill ) and therefore more gener- 
ally there are 0.6tot -1 mergers between z = 2.23 and z = 0.4, 
corresponding to O.Itqt -1 mergers per Gyr, where To = 1 Gyr. 



5 DISCUSSION & CONCLUSIONS 

The HiZELS narrow-band Ha survey selects star-forming galaxies 
within four well-defined volumes at z ~ 0.4 — 2.2 and flux limits 
with an SFR indicator which is unbiased in terms of stellar mass 
and is independent of its determination. In this paper we have used 
these properties to understand the star-forming population and its 
merger rate to help illuminate the processes responsible for the up- 
turn in the SFRD with redshift. 
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Figure 8. Left: Fraction of M20 identified mergers versus stellar mass for the four HiZELS redshift slices. Centre: Fraction of M20 identified mergers versus 
SFR for the four HiZELS redshift slices. From these plots we can see that the merger fraction depends on mass and perhaps SFR with the most massive and 
most star-forming galaxies having the lowest merger fractions. Right: Fraction of M20 identified mergers versus sSFR for galaxies with ENSFR > 0.2 for 
the four HiZELS redshift slices. This suggests that major mergers can lead to galaxies having unusually high sSFR compared to the typical value at a given 
mass and redshift. 




Figure 9. Left: Merger rates for the HiZ ELS sample ab ove a given mass against redshift. Fo r comparison, we include merger rates derived from: cl ose pairs 
jLin et al.l20o5 Lin08Ll 1); Gini/M 20 JLotz et alkooi Lotz08Ll 1); and galaxy asymmetry, JConselice et al.l2003Ll2009l:lL6pez-Saniuan et al.l2009i labelled 
C03, C09L1 1 and LS09L1 1 respecti vely) . The L I 1 denotes that these merger rates were origi nally sourced fro m their respective papers but have been corrected 
to the timescales calculated bv lLotz etall (201 1 ) using the galaxy evolution models of lSomerville et al.l J2008h . The samples of Lin08Ll 1, Lotz08Ll 1, C09L1 1 
and LS09L11 are all at M* > 1O 1O M while C03 is M > 1O 9 M . Right: The merger rates for HiZELS galaxies with M* > 1O 1O M above a given epoch 
normalised star formation rate (ENSFR = SFR/SFR*(z)). The points are offset by Az for clarity. From these plots one can see that there is no evidence 
for a significant evolution in merger rate when both the mass and the ENSFR of the galaxies are accounted for in the selection. 



By defining the epoch-normalised star-formation rate 
(ENSFR = SFR/SFR*(z)) we account for the increase in the 
typical star-formation rate of galaxies with redshift. In $2] we 
demonstrate that the number of galaxies above a given mass and 
ENSFR does not evolve significantly over the 6 Gyr from z — 0.4 
to 2.23. We also note that the HiZELS sample has already been 
shown to accurately trace the increase of the SFRD with redshift 
and that there is no st rong evolution in th e normalisation of the Ha 
luminosity function ds obral et al .1 120 1 3 ) . Taken, in combination 
this means the increase in the SFRD with redshift is not due to an 
increase in the number of star-forming galaxies of a given mass 
but instead must result from an increase in the amount of star 
formation in these galaxies. This can be describe d as an increase 
in the average sSFR for star-forming galaxies (Rodig hiero et alj 



2010: lElbaz etalJ koiD without a significant increase in their 
number density. Also, we note that the SFR* (derived from Z/H a ) 
evolves i n the same way as the typical sSFR for star forming 
galaxies telbaz et aDl201ll) . which implies that the luminosity of 
the knee in the Ha luminosity function is evolving significantly 
more rapidly than the characteristic mass of the stellar mass 
function. 

The size-mass relation for galaxies is assessed in £13.11 In 
order to do this for a large sample we need to use wide-field 
ground-based imaging. Hence we confirm that we can reliably re- 
cover the galaxy size determined from the HST CANDELS imag- 
ing by deconvolving the affect of atmospheric seeing from the 
ground-based imaging. We find that the size-mass relation is sur- 
prisingly constant out to z = 2.23, in agreement with the findings 
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of iBarden et"al1 (l2005h:llchikawa et all d2Q12h and at odds with the 
results of lTruiillo etaD d2007h : iMosleh et all d201ll) . The lack of 



strong size evolution at a given mass and the universal size-mass 
relation for star-forming galaxies in the range 0.4 < z < 2.23 sug- 
gests that this population have not experienced significant size evo- 
lution, through mergers or star formation, during this period. Any 
evolution that does occur must thus act to move the galaxy along 
the locus of the relation. The slope of this relation is also shallow 
and thus low mass galaxies are not dramatically smaller than their 
higher mass counterparts. Even if there is no direct evolutionary 
connection between the galaxy populations at each epoch then this 
lack of change in typical size suggests a universal evolution sce- 
nario. 

In order to study the merger rates of the HiZELS galaxies 
we test the Gini and M20 coefficients. By investigating these au- 
tomated methods of determining merger classifications we find 
that s EXTRACTOR parameters that define the segmentation map 
employed in these analyses are the most important factor in how 
well the method performs (see Appendix [A). We find that, for 
the segmentation maps generated by our set of S EXTRACTOR pa- 
rameters, the best delineation between mergers and non-mergers is 
M20 = —1.45 while the Gini coefficient provides no useful infor- 
mation. We acknowledge that other authors have found this not to 
be the case with M20 an d Gini being equally important in mor- 
phological classification (lLotz et alJ [20041 : IWang et al.ll2012h but 
we assume this is due to the differences in the construction of the 
segmentation maps (see Appendix lAt and potentially minor varia- 
tions in the normalisation of the M20 and Gini values, depending 
upon the exact nature of the morphological code used. The M20 
coefficient is found to be sensitive to mergers d o wn to a mass ra- 
tio of rsj 1 : 10 (in agreement with lLotz et al.l (l2010h ). We note 
here that not all mergers are star forming and as such we will 
miss 'dry' mergers which do not induce activity, although obvi- 
ously these will not be major contributors to the SFRD. As with the 
sizes we find that it is possible to use this morphological classifi- 
cation on ground-based data affected by atmospheric seeing, after 
applying a calibration derived from galaxies that are observed with 
both ground-based telescopes and HST. 

For the sample as a whole, without accounting for the Ha flux 
(SFR) limit or the increase in the sSFR of the star forming galaxies 
with redshift, we find that the merger fraction anti-correlates with 
both stellar mass and SFR. By combining these two results we find 
that the merger fraction correlates strongly with sSFR. This sug- 
gests that the more rapid the star formation is, the more likely it 
is to be driven by violent major mergers than secular processes. In 
fact we find that, ~ 50%, of starburst galaxies in our z ~ 2 sample 
have major merger morphologies. Therefore to achieve such high 
sSFR, these galaxies are undergoing major merger driven and not 
'main- sequence' star formation. Interestingly we see no evolution 
in the merger fraction of starbursts with a constant ~ 40 — 50% 
across all redshifts which suggests that merging is a universal pro- 
cess that can lead to a galaxy having enhanc ed sSFR for their epoch 
(iHopkins et a3200dlKartaltepe et all2012h . 

Finally we consider the merger rates of star-forming galaxies 
initially only limiting our sample on stellar mass, where some pre- 
vious studies have seen the characteristic merger rate increase to 
z rsj 1. However these other studies use photometrically selected 
samples where the method of determining stellar mass is directly 
linked to the determination of the SFR. As these two parameters 
are independent in HiZELS we can also select on SFR for a fair 
comparison across the redshift range, while also accounting for the 
increase in sSFR for typical star-forming galaxies with redshift. By 



applying these selections we see little evidence for an increase in 
the merger rates of typical galaxies over the redshift range consid- 
ered. Therefore even though there is an order of magnitude increase 
in typical SFR across the redshift range of our study this is not re- 
flected in the merger rate. This is strong evidence that it is not ma- 
jor mergers that drive th e increase in the SFRD w ith redshift, in 
contrast to the models of lSomerville etallfc oO 1 ) orlHopkins et alJ 
( 2006 1) and as o bserved in part bv lConselice et al.l (I2003L l2008h and 
iLin et alJ (2008) who find s ome evidenc e for an increase in merg- 
ing. Our result agrees with ISobral et all (120091) who find that the 
increase in SFRD between z — and z — 0.84 was primarily due 
to regular (non-merging) galaxies. 

Depending on the timescale r for which it is possible to view 
a galaxy undergoing a major merger using the M20 parameter we 
find that star forming galaxies with mass > 1O 1O M undergo 

- O.6T0T- 1 (3 if r = 0.2 Gyr) mergers between z — 2.23 and 
z — 0.4, corresponding to ~ O.Itot -1 (0.5 if r = 0.2 Gyr) 
mergers per galaxy per Gyr, where to = 1 Gyr. From analysis 
of the mass function of galaxies in COSMOS at z = 0.35-0.75, 
IPozzetti et alJ |2010) find merger rates of ~ 0.1—0.4 per galaxy per 
Gyr for galaxies with masses ^ 10 10-5 — 1O 11 M , in reasonable 
agreement with our findings (both are very sensitive to the c hoice 
of r). From a theoretical point of view Hopki mTet al.l d2010h com- 
pile data from a number of simulations and models (see references 
therein). The predicted number of mergers per galaxy with mass 

- 10 10 - 10 n M o per Gyr is found to increase with redshift from 
a value of ~ 0.05 at z = 0.4 to ~ 0.25 at z = 2.2 apparently 
lower than the values we find. Again, this is dependent on r so we 
are unable to provide solid constraints. 

In summary we find that the increase in SFRD is due to an 
increase in the sSFR of typical star-forming galaxies. The process 
responsible for this increase is not major mergers as we find that the 
merger rate does not increase in step with the SFRD. We therefore 
conclude that secular processes such as disc instabilities and/or an 
increase in the effective fuel for star formation are the main driver 
of the incr ease in the SFRD with redshift as predicted or observed 
bv others (Keres et al .1120051: iDekel etaDl2009l: iBower et al.ll2006l : 
iForster Schreiber et al.l 1201 ll ICacciato et al.ll2012l) . Although we 
note that it could also be driven by an increase in the minor merger 
rate (mass ratios < 1 : 10) which this study is not sensitive to. 

We also find a constant merger fraction for starburst galaxies, 
in that around half are major mergers across all redshifts, demon- 
strating that extremely violent events are required for a galaxy to at- 
tain enhanced sSFR for their epoch and leave the 'main- sequence'. 
Bringing these results together along with the lack of size evolution 
since at least z = 2.23 we can say that many of the properties of 
star forming galaxies are surprisingly constant over the ~ 6 Gyr 
covered in this study. 
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APPENDIX A: M 20 SIMULATIONS 

We discuss here the effect of the SEXTRACTOR property DE- 
blend_mincont on the segmentation map and derived M20 
value. In £13.2.11 we find that by analysing the M 2 o values and 
visual classification of a sample of z — 1.4 — 2.5 star form- 
ing galaxies, a segmentation map generated with a value of 
deblend_mincont=0.1 provides a demarcation between major 
mergers and non-mergers, with a boundary found at M 2 o ^ — 1.5 
(we later fix this value to —1.45). 

To quantify how the M20 value relates to a major merger we 
run some very basic simulations. The simulations comprise of mov- 
ing one artificial galaxy towards another and plotting the variation 
of M20 with distance. The artificial observations are created us- 
ing the GALFIT software with both galaxies being face on discs 
(i.e. Sersic index, n — 1) of the same magnitude and half-light 
radius. We first perform an analysis appropriate to the high red- 
shift star forming galaxies for which DEBLEND_MINCONT=0.1 is 
found to efficiently select mergers. For the simulations we use the 
appropriate values of the sky noise, magnitude zero point and PSF 
of the observation we are simulating. The results of this simula- 
tion are presented in the black points on Figure IA1I This demon- 
strates that for a DEBLEND_MINCONT=0.1 and CANDELS HST 
data the M 2 o value remains low and consistent with being a non- 
merger for galaxy separations down to ^ 1.6 arcsec (^13 kpc at 
z = 0.84 — 2.23), as the galaxies have individual segmentation 
maps. Once this separation drops below this value the two galaxies 
share the same segmentation map and thus the M20 value jumps 
dramatically as the top 20% of the light is now spread over two 
locations rather than one. As the separation decreases further this 
value lowers until a point is reached where the top 20% of the light 
is essentially co-located at the centre of a single bright galaxy and 
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Figure Al. The M20 value plotted against separation derived from a simple 
simulation of two identical face-on, disc galaxies approaching each other, 
as described in the text. The simulated ground based data are represented 
by red squares and the simulated HST data are black circles. To make 
the separations at which the M20 value jumps to be the same, we adopt 
a deblend_mincont=0. 10 and 0.03 for the creation of the space- and 
ground-based SEXTRACTOR segmentation maps respectively. The horizon- 
tal line at M20 = —1.45 represents the boundary between mergers above 
and non-merger below which we adopt throughout the paper. 

as such the M20 curve resembles a 'shark fin'. The distance over 
which this system would be classed as a merger is then ^ 1 arcsec 
which at z = 0.84 — 2.23 corresponds to a distance of ~ 8 kpc. 

Potentially the most important factor influencing the M20 
value for a given galaxy is whether it is derived from the space- 
based HST data as discussed above or from the ground-based imag- 
ing with a significantly larger PSF and different background charac- 
teristics. We first investigate this by using the original space-based 
value of DEBLEND_MINCONT=0.1 on the ground-based data and 
find that this gives a factor of ~ 2 larger range in separation over 
which the galaxies would be classified as a merger and would thus 
result in an increase in merger numbers relative to the HST imaging. 
A value of DEBLEND_MINCONT=0.03 accounts for this difference, 
equalling the separation over which a 'merger' occurs with the re- 
sults plotted as red squares in Figure IA1I This plot confirms the 
slope in the relation seen between ground and space-based derived 
M20 seen in Figure[5]and used to calibrate the ground-based values. 

For the lower redshift z = 0.4 sample, this angular distance 
range corresponds to 5.3 kpc and as such may miss some galaxies 
that would have been classed as mergers in the higher z samples. 
We therefore alter the value of deblend_mincont to 0. 1 1 for the 
CANDELS and 0.04 for ground-based imaging to account for this 
so that the same separation in kpc is used at each redshift. 

By varying the relative magnitudes of the galaxies and assum- 
ing the flux is linearly proportional to the mass and the size is pro- 
portional to the square root of the mass we test what mass ratio of 
mergers can be seen with this method. The result is that the M20 co- 
efficient is sensitive to mergers with a luminosity (mass) ratio dow n 
to ~ 1 : 10 (in agreement with the simulations of lLotz et al]|2010h . 
For mass ratios less than this the M20 coefficient does not increase 
significantly when the two galaxies share the same segmentation 
map. 



